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(ChPT) low-energy constants /la, /t. We found that in the instanton vacuum model the constant l-j 
is very sensitive to the shape of the instanton and the instanton vacuum parameters. We evaluated 
the constant l-j for two different zero-mode profiles and as a function of the average instanton size 
p and inter-instanton distance R. Our result agrees with an old "order of magnitude" estimate of 
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mass difference is small, ^ 1% of the observed experimental value. 
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I. INTRODUCTION 



The spontaneous breaking of chiral symmetry (S^SB) is one of the most important phenomena of hadron physics. 
It defines the properties of all the light mesons and baryons. Using the general idea of chiral symmetry, it was 
proposed in [l[ to use a phcnomcnological lagrangian, which has a form of the infinite series in the pion momenta 
and mass M^. The low-energy constants of the series expansion (LEC's) are the free parameters which encode 
the low-energy physics in a mo del- independent way. Up to now they were extracted phenomenologically from the 
experimental data, or from the lattice calculations ((MILC, ETM, JLQCD, RBC/UKQCD, PACS-CS)!!-!! within 
so-called Chiral Perturbation Theory (ChPT). 

One of the low-energy constants Ij is particularly interesting since it encodes the "pure QCD" part of the SU (2) 
isospin symmetry breaking (i.e. part which is due to u— and d— quark current mass difference, m„ — to^)- For 
example, the QCD part of the pion mass difference m'^+ — rii^a has a form [l[ 

{ml+ -ml,^)^^^^ —h{mu-mdf , (1) 

where B and F are the leading order parameters in the chiral lagrangian, and mu,md are the current quark masses. 
While experimentally the isospin breaking effects are known to a very high precision, separation of these effects on 
the "pure QCD" and electromagnetic parts has ambiguities and has been a subject of intensive debates From 
phenomenology the constant is known only with an "order of magnitude" estimate [l| , 

^7 - 5 X 10"^ (2) 

For this reason it makes sense to estimate this contribution in the framework of a reliable model. 

QCD instanton vacuum model, often refcred to as the instanton liquid model, provides a very natural nonpertur- 
bative explanation of the S^SB jlOl421| . It provides a consistent framework for description of the pious and thus may 
be used for evaluation of the low energy constants. Due to instanton-induced nonlinear interaction all the quark and 
meson loop integrals are regularized by the natural scale fj. ^ ^ 600 MeV in Pauli-Villars scheme [l^, where p is 
the average size of the instanton. This means that all the scale-dependent quantities, such as the quark condensate 
(qq) = (uu) + (^dd'^ and the difference 5 (qq) = (uu) — (dd), are given at the scale fi. Remarkably, the constant does 
not depend on the scale fj,. Recently [g^] it has been shown that this approach is able to give results consistent with 
phenomenological and lattice estimates for the constants I3J4, providing current quark mass dependencies of the pion 
mass TOjr and pion decay constant F^. 
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In this paper we would like to apply the instanton vacuum model for the evaluation of the constant l^. We extract 
the constant from the correlator {P^{x)P^{0)) using the relation [l[ 

f d'xe^'^^P%x)P°m = 4^ + ^ il") = ^-^^pi^^l, + O {m,q') . (3) 
J — q 9 

From the Eqn. ^ we may see that evaluations may be done in the limit m = ™"+™'' — > 0, and make only expansion 
over 

Sm = {niu - nid) ■ (4) 

The paper is organized as follows. In Section IIIII we discuss the general framework used for evaluation and write 
out the next-to-leading order (NLO) gap equation in the presence of the current mass split dm, which are needed 
for evaluation of the dynamical mass split 5M = M„ — Md- In Section ITVl we write out explicit expressions for the 
quark and meson propagators. In Section |V] we evaluate the effects of the mass split Sm on the quark condensate, 
S{qq) = {qq)u ~ {qq)d and extract the constant ^13. In Section IVll we evaluate the correlator {P'^{x)P^{0)) and extract 
the constant Ij. In Section IVlII we discuss obtained results, their uncertainty limits and draw conclusions. 



II. INSTANTON VACUUM MODEL 



The instanton vacuum model is based on the assumptions that the QCD vacuum may be considered as a dilute 
gas of instantons and antiinstantons, and the number of colors Nc is asymptotically large, Nc 00 (see the reviews 
[l5l[2^). While in general the sizes and local density of the instanton gas may be arbitrary, inter-instanton interaction 
stabilize these parameters. As it has been discussed in [2^ . the l/A'c-suppressed corrections due to the finite size 
distribution are indeed quite small, even for iVc = 3. Phenomenological, variational and lattice estimates lead to 
average instanton size p ^ 0.3 fm and inter-instanton distance R ^ 1 fm jl2l |. 

The partition function in the field of external scalar and pseudoscalar currents s = (sq -|- sr) and p ~ {po + pr) has 
a form [H 



Zn[sq,<7,sq,p,s,po] = y (iAexp(-re//[so, A,cr,so,p,s,po,i?i),??D,w]) , (5) 
T^ff = 5 + r™//, (6) 

S = ^ In A + 2 / d^xV - Tr In P + + " + + t ''/^ ' ^""^ , (7) 

V J ^ ' y p + isa + s-T + P075 + ip-T-f5 y ' ' 

The nonlocal formfactors F(j)) in the meson-quark interaction vertices come from the instanton-induced nonlocal 
interactions. Together with the factor [p + iso + s ■ f + P0J5 + ip ■ T75) in denominator, which subtracts the divergent 
high-frequency modes, they guarantee finite results for all the obscrvables in the instanton vacuum model. As it was 
discussed in [l^, [l^l , the divergent high-ferquency modes are responsible for renormalization of the parameters of the 
model. In what follows, we will fix them at the scale /i ^ p"^ ^ 600 MeV in the Pauli-Villars scheme [isj. 
The meson-loop correction y^ff to the effective action is given as 



F"/; [m, \,a]^-Tr\n( 4(5„ + ^Tr ( — ^ F^P) p 

"^f^ ' ' J 9 \ ^ cr2 \p + iso + S-T+po"f5 + ip-T"f5+icF<i^-TF ^' 



p + iso + s ■ f + po75 + ip - + i c F(f> • FF ■' 



$ • F = (^(7 + i75f'0 + if(T + 75?7j , (9) 

where c(A) = -^^^^^-^,5^ = ''^2N^^i^'' ^ color factor, F = {1, 75, ir, 17^75} is a set of matrices corresponding to 
quantum numbers of mesons present in the model, and we will use for the corresponding components of the field $ 
the notations $ = {<t, ?y, a, 0}. In contrast to NJL model, the variable A is a dynamical degree of freedom but not the 
parameter of the lagrangian. The current masses of the quarks come into play via constant external currents, viz. 

TO,, 4- md 

so = (10) 

TO,, — TOrf 

S3 = " • (11) 
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Notice that with respect to chiral transformations, the mesons may be separated onto two independent doublets 

^(7, and (77, (?) . The first doublet (^a, corresponds to the pion field U ~ (uq, u) in the notations of and the 

second doublet (77, a) is an additional degree of freedom which is absent in the chiral lagrangian. Now we are going to 
demonstrate explicitly on the example of the constant Ir that this additional degree of freedom (r/, a) gives an essential 
contribution to the constant l^ . As usual, the external currents (sq, s, Pq, p) generate nonzero vacuum averages of 

the fields (a) = t?^, (ry) = rjy and (^a, (j)^'^ ~ U = {uq, u). 

Due to the chiral symmetry expansion of the T^f / yields the general structure 

Teff[X,<J,s,pQ,ay,rjy,Ui] = Teff[m,X,a,s = 0,po = 0,ct„ = 0,riy = 0,Ui = 0] (12) 

+A (^{duof + {duf^ + B {sqUq +pu) +C (sqPo + ps)^ + (sor/^, + pa^f + a (pi + s^) + b {por]y + say) 

+c (ril + ay) + d {uopo + usf + e (mqPo + us) {u^-qy + uay) + / [m-qy + uayf + O {s^,p'^) , 

where we omitted the terms containing derivatives of the fields, since the external currents are constants, the constants 
A — V^ a — f should be evaluated with account of NLO corrections. The vacuum equations which follow from (fT2|) 
are 



dreff[m,X,ay,s = 0,po = 0,ay = 0,r]y = 0] _ dTeff[m,X,ay,s = 0,po = 0,ay =0,r]y = 0] _ 

dX day 
dVeff [m,X,a,s,po,ay,r]y, Uj] _ dVeff [m, X,a, s,po, ay , riy , u^] _ dVeff [m, A, g, s, po , ay , rjy , Uj 



day., 



drjy 



di 



(13) 
■ = (14) 



The coeSicients A, B are relevant for the 2-point correlators with intermediate pious and A ^ and B oc M"^ in 
•niu = TTLd limit. The constants B, C, V are irrelevant to our problem since they are constants in front of the term with 
chiral doublet x — iso,p) which we put to zero in the current paper. 

The Eqns (|13p are responsible for the dynamical mass generation and will be discussed in the next section. The 
Eqns (jl4p may be explicitly written as 



dTeff [m,X, a, s, po,ay,r]y, Uj 
drjy 

dTeff [to. A, cr, s, po , ay , 7?t, , Uj 
day^i 



= 6po + 2c-qy + euo (itoPo + us) + 2/uo (ito?7i, + uay) = 
= bsy + 2cay + eu (mqPo + us) + 2fu (u^rfy + uay) = 



Multiplying Eqn (jlSp on uq and Eqn. (1161) on u and adding results, we may find: 

h + e 



UQTjy + uay 



Repeating the same trick with po and s, we may get 



2(c + /) 



1 

Poriy + say = - — 
2c 



(uoPo + us) . 



.b + e 



and repeating the same trick with ?/„ and (?„, we may get 



2 I -'2 
Vy + <^y 



1 



2c 



6 + e 



2/ 



Combining results (|17m9p . we may get for the effective action 

b 



reff[X,a,s,pn,u,] = ... + « (po + ^ ^ 



'2(c + /) V2(c + /) 



6 + e 



(wopo + us) 



b{pl + ^) + \ e~ /— -7 (woPo + us) 



4c 



b 

4^ 



b + e 



+d (uqPo + Ms)^ — e 



6 + e 

c+/y '^4(c+/) ' •' \2{c + .f) 

b + e 



(15) 
(16) 

(17) 

(18) 

(19) 

(20) 



b + e 
e— ^ + ./ 



(uoPo + us) 



2{c+f) 



[uoPo + us)^ + f 



2(c+/) 



(uqPo + us] 



^...+ (a-^)(p^ + .-) + (d-£(e-/^)-e 



4c 



c+fj 4(c+/) 



(uopo + us)'^ + O (x, X^x) 
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where we omitted the terms which are proportional to the chiral doublet x- The terms shown in ((20)) are explicitly 

chiral invariant and correspond to the terms (x^x) sind (x^C/) in the chiral lagrangian Respectively, for the 
constant we may deduce 

= ^-V'- ^^\^) (21) 

Thus we can see that in addition to the term d in numerator there are three other terms which correspond to 
contributions of additional mesons. As we will see from the following sections, these contributions have different signs 
and approximately the same order of magnitude as the term d. The formula (|21[) proves that we have to consider 
correlators instead of direct comparison of the terms in the expansion of the lagrangian. 

Below we will not evaluate the constants A — V^a — f, but instead evaluate the correlators directly. 



III. GAP EQUATION 



The next-to- leading order (NLO) gap equations which follow from the effective action ^ have a form 



0-3 



A 



dS_ 
'da 
dS 
das 
dS 



dX 



where we used notations 



4^2- 




4^1- 


V 


N 

V ^ 


2V 


1 




+ 2^ 


J (27 



TrUSMip)T3S{p) 



;Tr [S{p) {iM{p) + ^T3<5M(p))^ 



y(gap)(ij] 


'(9) = 




f d'p 
(2^)4 








(2^)4 


fAgap){ij] 
^3 


'(9) = 




f d^p 

(2^)4 



Tr 



Tr 



Tr 



M{p) 



M{p + q) 



p + i/u(p) + iT3Sjj.{p) p + q + i^{p + q) + iT3Sfi{p + q) 



M{p) 



p + + iTsS^iip) 
Mip)SM{p)T3 



M{p + q) 



p + q + i^i{p + q)+iT36i-J.{p + q) 
M{p + q) 



(p + i^i(jp) + iT36fi{p)f I P + q + iKP + 9) + iT36fi{p + q) 



(22) 
(23) 

(24) 

(25) 
(26) 

(27) 



explicit expressions for the vertices (|25ll27p are given in Appendix [Xj and the propagators used for evaluation are 
written out in Section |IVl In general, equations (|22ll24p can be solved only numerically. 



A. Expansion over Sm 



For the special case when dm is small, it is possible to solve the equations (|23p making a systematic expansion over 
the small parameter 6m. For our purpose it suffices to keep just the first order corrections. From the first and the 
third gap equations in p3p and the structure of the vertices p5|26p we may conclude that the vacuum expectation 
values for (ct), (A) get corrections only in the second order over Sm, thus in the first order they remain the same as 
for Sm — 0. The equation for the (a^) has a form 



dS _ .2/ .2 « /• d^p {p'-fi^{p))M{p){Sm + eM{p)) 

o-3^«-4e {a) -8eNc / -77^- ;7-72 28) 

00-3 J i2TTp {p^ + ^J■^{p)) 

^ ^ '^%EK3^''^('?)%(9) = 0' (29) 



t2 



(27r) 



where e = 



5 



V^i^'^%(g) « ^ [ 4^MHp)M{p + q)x{ 



2^2 J (27r)4 





(n(.°j(g)-nil,(g)) 




-q) 


-S- 


ip)s- 


ip)s- 


b^ 


' l)\o{5m) + 




211(7(73 


Tr[S+{p)S+{p)S+{p^ 


-1) 






ip)s- 


b^ 


- l)]5m=0 ~ 


E 

1 = 1,2 


ni?(fc) 


Tr[S+{p)S+{p)S-{p^ 


-q) 


-S- 


ip)S- 


{p)s+ 


(p^ 


' 9)]o(5m) ~ 






Tr [S+ip)S+ip)S+ip- 


u) 


- S- 


-ip)s- 


-ip)S- 


[p- 


'^-']e'(c5m) ^ 




2n,^3(fc) 


Tr [S+{p)S+{p)S+{p- 




+ S- 


-ip)S- 


-ip)S- 


[p- 




E 


ni°^(fc) 


Tr [S+ip)S+ip)S.ip- 




-S- 


-ip)S- 


^ip)s+ip- 


^ ^-'l 0(5™) 



ij_ = l,2 

the superscript (0) on the propagators and subscripts on Tr [...]„ indicates that the proper propagator is to be taken 
in the hmit 5m = or just collecting the first C'((5m)-correction. One can notice that ([28| has a form 



e{Xe + y5m) = Q, (30) 

where the coefficients X,y should be evaluated with account of C'(l/A^c)-corrections. A trivial nonzero solution is 
e = —5my/X, which corresponds to 

- 6m (1 - M f{p)y/X) . (31) 
While in general case the explicit expression for the formfactor has a form 

Hp) - -a;^ {h{x)Ko{x) - h{x)K^{x))^^,_, , (32) 

in order to speed up the evaluations here and below we consider two parameterizations of the formfactor. The first 
one is a simple "dipolc" parameterization ([l^]) of the form 

/(P) = .\ 2 ' (33) 
2 + p^p'^ 

which coincides with p2p in the region of small p < 2/ p. The second parameterization has a form 

1 



f{p) 



(34) 



where the free parameters ai.-.as are fitted to ([32)) . and ([34]) agrees with ([32|) both for the small and asymptotically 
large p. We will refer to and ([M)) as dipole and quasibessel parameterizations respectively. Direct comparison of 
the two close parameterizations is important in order to demonstrate that the results of this paper are very sensitive 
to the shape of the instanton. 

We summarize results obtained for the constants X ^ y with different parameterizations of formfactor in Table H] 
As we can see, the 1/iVc-corrections X^lq and yNLO are small compared to X^o and yLO respectively, so the l/Nc- 
expansion works very well here. It is important to note that both in the leading order and next-to-leading order the 
mass Sp{p) changes sign for p ^ 0.5 GeV, so we have a compensation of the small-p and large-p regions. 



IV. PROPAGATORS 



In this section we would like to discuss the propagators of the quarks and mesons in presence of the mass split Sm. 
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Xlo X 10^ 


yLO X 10^ 


-{y/x)Lo 


Xnlo X 10^ 


yNLO X 10^ 


— (y / X)lo+nlo 


Dipole 


-8.50 


-3.53 


-4.16 


0.80 


0.61 


-3.80 


QuasiBessel 


-9.04 


-3.33 


-3.68 


1.02 


0.69 


-3.30 



Table I: Parameters of mass split in different parameterizations of formfactor. Dipole corresponds to (|33p . QuasiBessel 
corresponds to (|34|) . Dimensions: [Xlo.,nlo\ = [GeV\ [yLO,NLo\ = [GeV'% 



A. Quark propagator 



Here we consider only the leading order quark propagator, the NLO corrections to the quark propagator will be 
considered as separate meson loop corrections to proper correlators. Since the operator j3 + i/i(]3) + i5/i(p)r3 is diagonal 
in the flavour space, its inversion is quite straightforward, with 



S±{p) = , (36) 



where fJ.±{p) = /i(p) ± S^{p). 



B. Meson propagator 



For evaluations in this paper we have to evaluate the meson propagator with account of l/7Vc-corrections. However it 
is important to note that NLO evaluations are needed only for Hjj,j(0), Res Il^^{q) Jl,i^{0) ^ all the other components 

and expressions for g 7^ may be evaluated in leading order, which significantly simplifies the task. Due to Sm ^ the 
propagator is nondiagonal in indices (ij)-wc get additional transitions a a and i] ^ 0. Inversion of the propagator 
is trivial and gives: 

rr ^ )33 tt ^ )oO (o-j) 

(n-i)oo(n-i)33-(n-i)o3(n-i)3o' i^-')ooi^-')33- i^-')o3i^-'ho' ^ ' 

TT V /30 TT V /03 /oo\ 

°' (n-i)„„(n-i)33-(n-i)„3(n-i)3o' '° (h-i)„„ (h-i)33 - (h-i)„3 (h-i)3o' ^ ' 

where we used a shorthand notation (0,3) = (cr, a^) for positive parity mesons, and (0,3) = (77, (jy^) for negative 
parity mesons. 

1. Leading order 
In the leading order for the components (H~^)^^. we have 

(H-i)^^. - 4S,, + ^Tr (q{p)T,Q{p + q)T,) (40) 

where 

Q{p) = iM{p)S{p) ^ ;f_[^.^ , (41) 

p + ^^l{p) + idfj.{p)T3 

and explicit expressions for the components are given in Appendix El 
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2. NLO correction 

As it was discussed earlier, we need a few values for propagators in the next-to-leading order. Since the NLO 
evaluations arc numerically slow, from the very beginning wc will concentrate on evaluation of the following quantities: 



lim 

(ij->0.m-i-0,(5rn->0) 



nJ(g), 



lim 



(g-i.O,rn->0.i5m-i.O) 5m q^=0 



All the terms which do not contribute to one of these limits will be omitted. For the sake of brevity, below we use 
notation 



2 p + in(p) 



(42) 



For the pion propagator Il^,f,(q), we may use the chiral limit and put m,5m to zero. The NLO expression for the 
pion propagator has a form 



n 



(af))-l 



(9) 



1 f d^k 
1 



ny(fc) (2rrp (Q(p)i75r"Q(p + q)T,Q{p + q + k)T,Q{p + 9)1757 



~Trp {Q{p)i^5T''Qip + q)r,Q{p + q + k)ij5T''Qip + fc)r,)) 

4 r d^k 



(27r 

Trp {Qip)ij5r''Q{p - q)r,Q{p - q - k)T,) . 
In complete analogy, for the 77-meson propagator n,,,,(0) we have 



jn,{k)nj{k + q)Trp (0b)^75r"Q(p + q)T^Q{p + q + k)Tj) 



n-^(o) 



45-6 + -^Trp {Q{ph,Q{ph5) 



1 



%(fc) {2Trp {Q{p)-i^Q{p)T,Q{p + k)T,Q{p)^^) - 



(T^ J (27r)" 

+ Trp [Q{p)i^Q{p)T,Q{p + A;)75Q(P + k)T,)) 
4 f d^k 



n,ik)nj{k)Trp {Q{p)-i^Q{p)V,Q{p + k)Tj) X 



(27r) 

rrp(Q(p)75Q(p)r.Q(p-fc)r,), 

and again we can make evaluations in the chiral limit. 

The nondiagonal matrix element n^0(O) is O (Sm), so we will extract explicitly Sm and after that the evaluation 
of the constant will be done in the chiral limit. Evaluation is quite tedious since LO propagators have nondiagonal 
components. The corresponding expression has a form 



n-;(o)-n 



l]cj> 



l]cj> 



(0), 



where 



(43) 



T^Trp {Q+{p)ij5Q+{ph5) - T^Trp {Q~{p)i-fzQ^{p)-fz) 



2cr2 



(44) 



n 



(1- 

ri4i 



s)-l 



(0) 



d'^q 



{l — mes. 
'^3 



d^q 



d'^p 



(27r)4 ^ *J J (27r)4^ '''^'2cr^J (2tt)^ 

ij ij 

X Tr i2Qip)r^Qip)T^Qip)T,Qip + q)T, + Q{p)T,jQip)T,Q{p + g)r^Q(p + q)r,) , (45) 
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LO 


Mass Shift 


Mass Split 


Meson 


All NLO 


Total 


QuasiBessel 














n-,;(o) 


5.65 X 10"^ 


9.16 X 10"^ 





-8.88 X 10"^ 


2.86 X 10"^ 


5.93 X 10"^ 


-*n;;(o) 


-4.73 X 10"^ 


0.88 X lO^-' 


0.19 X 10"^ 


-2.30 X 10"=* 


-1.23 X 10"-' 


-1.89 X 10~^ 




1.24 X 10"^ 


-0.59 X 10^2 





0.11 X 10"^ 


-0.49 X 10"^ 


0.76 X 10"^ 


m 


2.03 X 10"^ 


-0.77 X 10"^ 





0.31 X 10"^ 


-0.45 X 10"^ 


1.58 X 10"^ 


B 


1.64 


— 





— 


— 


2.09 


Dipole 














n-,}(0) 


5.65 X 10^^ 


9.16 X 10"^ 





-9.72 X 10"^ 


-0.56 X 10"^ 


5.09 X 10^^ 




-1.72 X 10"^ 


1.47 X 10"^ 


0.17 X 10"^ 


-4.95 X 10"^ 


-3.32 X 10"^ 


-5.20 X 10"^ 


F'^ 


1.36 X 10"^ 


-0.54 X 10"^ 





0.31 X 10"^ 


-0.23 X 10"^ 


1.12 X 10"^ 


im) 


2.18 X 10"^ 


-0.72 X 10"^ 





0.31 X 10"^ 


-0.41 X 10"^ 


1.76 X 10"^ 


B 


1.60 











1.57 



Table II: In this table we give the numbers obtained for propagators and other relevant constants. is used for evaluation of 
Res (n^) , (gg) is used for extraction of constant B. 



n 



(2-mes)-l^n~l — 
ri<t> 



(46) 



y^k\<D = [Trp {Q{p)l^Q{p)T^Q{p + q)Vk)] , 
V}k \<l) = [Trp {Qip)il5r''Qip)T,Q{p + q)Tk)] , 



and explicit cxpessions for the verices contributing to (j45ll46p are given in Appendix [X] 

Numerial results of evaluation are presented in Table [III As we can see, even in the leading order (LO) there is 
a strong sensitivity of the propagator n~_^(0) to the shape of the instanton (formfactor F{p)). This dependence is 
discussed in more detail in Section IVlIl 



V. QUARK CONDENSATE 



Due to the mass split 6m there is a flavour difference for the quark condensate S (qq) = (uu) — (dd) . In the leading 
order this split is 



^ {<l<l)i 



-Tr {t^S{p)) = 4A^, 



d'^P ( M+(p) f^-{p) 

w.-v/^;;---cy ^^^^^ ^ + ^2^^p) p2 + ^2^^p) 

In the NLO evaluation is also quite straightforward, with 

^(99> meson = 



d^q 



V:f''\q)^- 



(2-) 
d^p 



(47) 



(48) 



-^M{p)M{p + q)Tr ^ !j5(p)r,5(p + q)r,Sip) 



for meson corrections plus corrections from mass shift and mass split {1/Nc corrections to Mq and M3), and explicit 
expression for p8)) is given in Appendix El 

Results of numerical evaluation are presented in the Table Hill As one can sec, due to the large NLO corrections to 
the mass split A£u(p) — M^lpYthe NLO corrections are larger than the LO result. 

Using formula (11.3) from [l[, it is possible to get for the constant an estimate [28| : 



{{uu) - (dd)) 



\LO+NLO 



0.10 (5m 

AB'^Sm 



5.48 X 10~^. 



(49) 
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LO 


Mass Shift 


Mass Split 


Mesons 


All NLO 


LO+NLO 


Dipole 


-0.20 


6.07 X 10"^ 


1.26 X 10~^ 


1.09 X 10~^ 


7.45 X 10~2 


-0.13 


QuasiBessel 


-0.18 


6.11 X 10"^ 


1.29 X 10"^ 


2.99 X 10"^ 


7.70 X 10"^ 


-0.10 



Table III: Different contributions to — LO: Leading order result Mass Split: Contribution due to NLO correction 

to mass split Mu(p) — Md{p)- Mesons: Contribution of mesons, All NLO: sum of contributions of mesons and mass shift, 
LO-fNLO-flnal result. 

Figure 1: Contribution to the (P^P")-correlator in the leading order. 




VI. EVALUATION OF THE CONSTANT h 

According to it is possible to evaluate the constant Ij from the correlator {P^{x)P^{0)) as 

P,iq) = [ d^xe^^^P^x)P"m = 4^ + O {q') = + O {m.q^) , (50) 

J mil:. — q'^ q'^ — 

where the constant B is one of the phenomenological parameters of the chiral lagrangian (see Table |ll]) , the mass of 
the pion = in the limit m — > and m„, are the current quark masses. Since we are interested only in the 
residue of the correlator, we should consider only 1-particle reducible diagrams with pion in the intermediate state. 

In the leading order, there are two diagrams shown in the Figure [TJ Obviously, only the diagram on the right-hand 
side contributes to the residue, yielding 



^a^^l?) = E if°('Z)^'i^°(9)n„-(q) = Lr,{q)L^{q) (n„„(g) + H^^q)) + (L^/g) + Ll{q))IV^^{q) , (51) 



where 



i^^lg) = -iy^zM/(p)/(p + q)[rr(5+(p)755+(p + g)75))+Tr(5_(p)755-(p + g)75))]= (52) 
= 4z7Vc/ -0^,Mf{p)f{p + q) 

Li^'iq) = II ^Mfip)fip + q)[TriS+{ph5S+ip + qh,))-TriS^{p)^,S^{p + q)^,))]= (53) 

+ p-q + ^J.+ {p)^^+{p + q) p^ + p-q + {p)fj.- {p + q) 

(p2 + ^Kp)) ({p + q)2 + ^l(p + g)) (p2 + ^^(p)) ((p + qy + + q))y 

and we used identities 



p'^ + p ■ q + ^J.+ [p)p+{p + q) ^ p'^ + p ■ q + n- {p)fi- {p + q) 
(p2 + (p)) ((p + g)2 + (p + q)) (p2 + (p)) ((p + q)2 + (p ^ 



= - 4iV,y" -l^Mf{p)f{p + q) 



i?^0(g) = Lf{q), (54) 
i?^0(g) = L^O(9). (55) 

In the next-to- leading order there are seven diagrams shown in the Figure [21 Obviously, only the diagrams 3-6 from 
the second and the third row contribute to the residue in pion pole. The explicit expressions for the corresponding 
diagrams are given in Appendix]^ Using p7l38p . one may immediately get 



ResU^^q) w -U-^{0)n^{0)Resn^^{q),+O {Sm\ m) (56) 
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Figure 2: Contribution to the (P'^P'')-correIator in the next-to-leading order. 





LO 


Mass shift 


Mass split 


Meson 


All NLO 


Total 


QuasiBessel 














— zZ/77 


4.07 X 10"^ 


-1.53 X 10"^ 





-2.33 X 10"^ 


-3.87 X 10"^ 


1.94 X 10"^ 




-0.15 X 10-=* 


-8.03 X 10"^ 


-2.34 X 10"^ 


-6.93 X 10"^ 


-1.73 X 10"^ 


-1.74 X 10"2 


h 


0.17 X 10"** 










1.198 X 10-* 


Dipole 














— iLr, 


4.35 X 10"^ 


-1.45 X 10"^ 





-2.38 X 10-2 


-3.83 X 10-2 


5.22 X 10-^ 




6.71 X 10"^ 


-1.06 X 10"^ 


-2.10 X 10"^ 


-1.70 X 10"^ 


-2.97 X 10-2 


-2.30 X 10-2 


h 


0.34 X 10"^ 










1.00 X 10-^ 



Table IV: Evaluation of the residue Res^2^_^2 {P3P0). See Eq. (|57p for more details on meaning of Ln,L^. The first column 
is the LO result, columns 2-5 are NLO corrections, column 6 is the total result. In columns 7-8 we give results for h in LO and 
NLO (See the Table HH for numbers used in evaluation). 

So in evaluation of the residue Res {P3P0) ~ Res (P3P0) one has to keep only the terms 

q^— — m^ q^—Q 

{P3P0) = = L^{q)L^{q)I^4>M + {L^M) + L^(g))n^4g) + non - singulars, (57) 

all the other terms which are not written out explicitly do not contribute to the residue. 

Results of numerical evaluation are presented in Table ITVl As one can see, the model is extremely sensitive to the 
change of formfactor. The reasons of such strong dependence will be discussed in Section fVIIl 

VII. CONCLUSION 

In this paper we evaluated the effects of the current quark mass split on the dynamical mass, quark condensate 
and correlator (PsPq)- From these data we extracted the low energy constants hz,h- We found that the dynamical 
quark mass 5M is negative, so as one can see from the left pane of the Figure [3l the momentum-dependent mass 
(5/i(p) = 5m + 5Mf'^{p) has different signs for small and large momenta. Due to cancellation of these contributions, 
we got very strong sensitivity of all quantities discussed in this paper to the details of the instanton vacuum model, 
such as the shape of instanton (which comes via the formfactor) and instanton parameters. In the right pane of 
the Figure [3] we demonstrate explicitly this fast dependence on the example of the leading-order integrand of ^^(0). 
As it was explained above, due to different signs of large and small-p contributions, we have partial (solid line) or 
almost complete (dashed line) cancellation, which leads to the strong dependence on parameters of the model. Similar 
behaviour is observed for all quantities where the dynamical mass split 5M(p) contributes, both in the leading and 
in the ncxt-to-leading orders. 

It is necessary to note that the instanton vacuum model contains chiral doublet (77, {7)-additional degree of freedom 
which is absent in the chiral lagrangian, and the cancellation of the different contributions is due to the dynamics of 
the field. If we set —y/X = in (|3T|) and thus effectively eliminate the contribution of the 0-3, we can see that the 
dynamical mass split (5/x(p) is constant for all momenta p, and cancellation of different regions does not happen. 
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Figure 3: Left: dependence of the dynamical mass split Sfi{p) = 5m + 5Mf^{p) on the quark momentum p. Right: Instanton 
shape dependence of the integrand of L^{0) in the leading order. g{p) is the integrand of the Ean l53l 




Figure 4: [Color online] Dependence of the constant Ij on the instanton vacuum parameters p and R 



One of the consequences of the above-mentioned sensitivity of to model details is that uncertainty of the instanton 
vacuum parameters (average instanton size p and inter-instanton distance R) leads to increased uncertainty in the 
final prediction for Zy. As it has been discussed in (22l |. different methods estimate the model parameters are in the 
range p ~ 0.32 — 0.35 fm, R ^ 0.8 — 1 fm. While the uncertainty in p, R is just ^ 10% and is unimportant for most 
evaluations, for the constant it leads to sizeable uncertainty in the final result. Using the Figure SI we may get for 
Ir an estimate 



/y - (6.6±2.4) X lO""*. (58) 

The result (j58p agrees with a phenomenological estimate ^ within uncertainty limits. Using ([1]), we may obtain 
for the pure QCD contribution to the pion mass difference 



{mU-mlo) r^lAxlO-^GeV^, (59) 
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i.e. ^ 1% of the experimentally observed difference. This result does not contradict the well-known fact that the pion 
mass difference has electromagnetic origin j25l - l27j . 
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Appendix A: Explicit expressions for some vertices 

In this section for the sake of completeness we would like to present some explicit expressions for the meson-quark 
interaction vertices which are used in this paper. For the quark-meson vertices (|25ll27p which come into the gap 
Eqns. (|22ll24p we may get 



2 '''■^> a^J (27r)4 \p+iti{p) + iTs5fiip) 



1 f d^P M{p) 

M{p + q) p. 

p + q + ifi{p + q) + 2r3(5/i(p + q) ■' 
^ ^ '^'P-M{p)M{p + q)x{ 



2a2 J (27r)4 

(n,,(g) - n^.^M) Tr [S+{p)S+ip + g) + S-{p)S-{p + q)] + 
2n,,3 Tr [S+ip)S+{p + q)- S^{p)S^{p + q)] - 

n,.(fc) Tr[S+ip)S-{p + q) + S-{p)S+ip + q)]- 



= 1,2 



{U^^ik) ~ n^3^3(fc)) Tr [S+{p)S+{p + q) + S-{p)S-{p + q)] 
2n,^3 (k) Tr [S+{p)S+{p + q)-S^ {p)S^ ip + q)] 

n^.(^) Tr[S+{p)S.{p + q) + S.{p)S+{p + q)] 

±=1,2 

} 



M{p + q) ^ 

p + q + i^ip + q) + iT^6i^i{p + q) ^ 

^ ^ '^'PrM\p)M{p + q)x{ 



2a^ J {2tt)^' 

{n,,{q) - Il^,aM) Tr [S+{p)S+{p)S+{p + q) + ^_ (p + q)] + 
2n,,3 Tr[S+{p)S+{p)S+{p + q) - S-{p)S-(p)S-{p + q)] - 

Y ^a.{k) Tr [S+{p)S+{p)S^{p + q) + S^{p)S^{p)S+{p + q)] - 

■ ±=1,2 

(n,,(fc)-n^3^3(fc)) Tr[S+{p)S+{p)S+{p + q) + S-{p)S-{p)S-{p + q)] - 

2n,^3(fc) Tr [S+{p)S+ip)S+ip + q) - S-ip)S-{p)S-ip + q)] + 

J2 n^.(fc) Tr[S+{p)S+{p)S^{p + q) + S.{p)S^{p)S+{p + q)] 

■±=1,2 

} 
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p + g + + g) + iT3Sfi{p + q) ^ 

i / (^^^ b)^^^(P)^^ C?' + 9) X { 
(n,,(<7) -n<,3,3(q)) Tr[S+{p)S+{p)S+{p + q) - S-{p)S-{p)S-{p + q)] + 
2U,,, Tr [S+{p)S+{p)S+{p + q) + S-{p)S-{p)S-{p + q)] - 

J2 n,.(fc) Tr [S+{p)S+{p)S-{p + q)- S-{p)S-{p)S+{p + q)] - 

(n,,(fc)-n^3^3(fc)) Tr[S+{p)S+{p)S+{p + q)-S^{p)S.{p)S^{p + q)] - 
2U,^,{k) Tr [S+{p)S+{p)S+{p + q) + S^{p)S^{p)S^{p + q)] + 

n^.(fc) Tr[S+{p)S+{p)S-{p + q)-S-{p)S-{p)S+{p + q)] 

i±=l,2 

}■ 

For the components of the leading order meson propagator (|40p we may get the foUowing expUcit expressions (29j 
(n-i)^^ = i+^Tr{Q+{p)Q+{p + q)) + ^Tr{Q^{p)Q^{p + q)) (Al) 
(n"').3.3 = 'i-^Tr{Q+ip)Q+{p + q))-^Tr{Q^{p)Q^{p + q)) 
(n"')<.„.3 = ^Tr{Q+{p)Q+{p + q))-^Tr{Q^{p)Q^{p + q)) 

= 4-^rr(Q+(p)g+(p + g)) -^rr(g_(p)g_(p + g) 
(n"')0303 = 4 + ^rr(Q+(p)g+b + g))+^rr(g_b)g_(p + g) 

(n-^)^^3 = -^Tr[Q+{p)Q+ip + q))+^Tr[Q^{p)Q^{p + q) 

= 4 + ^ ^Tr(g„(p)r,g„(p + g)ro, j) ^ (o, 3) 

where Q±ip) = ^^^^^^^^^^(p) , Q±{p) = ~1^Q±{p)ib = p-^,l!^)%^(p) ■ 
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The 1-loop correction to the propagator n|^^ ^ (0) has a form 



ij 



J j^Tri2Qip)T,,Qip)T^Qip)T,Q{p + q)T,+Q{p)r^Q{p)T,Q{p + q)r^Qip + q)T, 



4 

i4 



(27r 



* ^ ^-M\p)M{p + q)x{ 



2cr4 7 (27r)4' 

(nW(g) - U^ZAi)) Tr (-2M(p) (5+(p)5+(p)5+(p)5+(p + q) - S^{p)S^{p)SAp)S^{p + g)) + 
+M{p + q) {S+{p)S+{p)S+{p + q)S+{p + q)- S-{p)S-(p)S-{p + q)S-(p + q)))^^^^) + 
2U,„,{q)Tr{-2M{p) {S+{p)S+{p)S+{p)S+{p + q) + S-{p)S-{p)S-{p)S-{p + q)) + 
+M{p + q) {S+ip)S+{p)S+{p + q)S+ip + q) + 5_ (p + q)S-ip + <z))),,„^o - 

^ nW(^)rr (-2A/(p) {S+{p)S+{p)S+{p)S-{p + q) - S-{p)S^{p)S-{p)S+{p + q)) + 

i±=l,2 

+ M{p + q) {S+{p)S+{p)S-{p + q)S-{p + q)- 5_ + q)S+{p + + 

(n(0)(g) - n^'^^lq)) Tr {2M{p) {S+{p)S+{p)S+{p)S+{p + q) - S^{p)S^{p)S^{p)S^{p + q)) + 
M{p + q) {S+{p)S+{p)S+{p + q)S+{p + q)- S-{p)S-{p)S^{p + q)S-{p + + 
2U,^,{q)Tr {2M{p) {S+{p)S+{p)S+{p)S+{p + q) + S-{p)S-{p)S-{p)S-{p + q)) + 
M{p + q) {S+{p)S+{p)S+{p + q)S+{p + q) + S^{p)S.{p)S^{p + q)S^{p + (?))), ,„^„ - 
ni°^(g)Tr {2M{p) {S+{p)S+{p)S+{p)S-{p + q) - S-{p)S-{p)S-{p)S+{p + q)) + 

M{p + q) {S+{p)S+{p)S-{p + q)S-{p + q)- S-{p)S-{p)S+{p + q)S+{p + q)))^^,^^ 
The two-loop correction to n,,0(O) has a form 

n!,;-"^^^"^(0) ^-±1 -0jn,(,)nH(g)F/;')(g)F;«(,), (A2) 

where 

Vi^\q) = [Trp {Q{p)-l^Q{p)T,Q{p + q)Tk)] , 
y^k\1) = [Trp {Q{p)^l^T^Q{p)T^Q{p + q)Tu)] • 
In exphcit form, with account of O ((5TO)-counting, (|A2| has a form 

Y^IV,,{qnui{q)V.^:^\q)V^f\q) = (A3) 

ijlk 

+ V^^l{q)V^ll{q)n,,{q)n^,M + Vjj^f^^{q)V^tl{q)n,,^M^,,{q) 

+ V^^ {q)V^t^ {q)Ii,, {q)IV,„ [q) + V^^^ {q)V^X {q)Tl^,^., {q)Ii^,^, [q) (A4) 

((7)n,,03(g) 

+ K.^,^^ (<?) ('7)n.. (g)n„^3 (g) + Fj,^) (q)F4^) (g)n..3 (g)n,, (g) (A5) 



where the vertices have an explicit form: 

= Trp {Q{p)l,Q{p)Q{p + g)*r375) = 



"^'^ M\p)M [p + q) Trp {S+ip)S+ip)S+ip + q) - S-{p)S-ip)S-{p + q)). 



2 J {2nf ^^^^-^^^ ' ■P\-^\i')-^yr)-^\t' ' h, ^-y^ i^-yt- i^-yt- . ^no[&^ 
y^t^^) = Trp {Q{p)il5T^Q{p)Q{p + g)»T375) = 

^ ^ '^^-M\p)M{p + q)TTj,{S+{p)S+{p)S+{p + q) + S-{p)S^{p)S-{p + q))^ 



y^^'^M) = (Q(p)75g(p)75g(p + 9)«r3) = 

-\\ -^^M\p)M {p + q) Trp {S+{p)S+{p)S+{p + q) - 5_ (p)5_ (p + q))^^,^ 

Vv^fsi<l) = {Q{p)il5T^Q{p)l^Q{p + q)iT^) = 



1 r d^p 



M\p)M {p + q) Trp {S+{p)S+{p)S+{p + q) + 5_ (p + 



2 J {2n f 

(7±0± 

1 /■ d'^p 



- , ---^M\p)M{p + q)Trp{S+{p)S+{p)S^ip + q) + S.ip)S.ip)S+ip + q))^^^^ 
^ J (27rj 

Vi"^'-. (g) Trp iQ{p)iT3^5Qip)iT±j5Q{p + q)iTj_^5) = 

' f '^^^-^M\p)M{p + q)Trp{S+{p)S+{p)S-{p + q)-S.{p)S.{p)S+{p + q))^ 



= Trp {Q{p)i^Q{p)Q{p + g)75) = -V^/^Sq) = O (-5™") , 
V^t^ (q) = Trp {Q{p)in^,Q{p)Q{p + g)75) = ^j;'^ (9) = 0{5v^), 

^^3^3 = (Q(p)*T375gb)*T3Q(p + g)*T375) - -V^tM) = ^ ('^"^) • 
Using the last four equations (|A15IIA18|) . the two- loop contribution (|A3[) may be cast into the form 

Y,n,M^,M)v^^\q)v}f\q) = 

ijlk 

v^;liq)v}^^^iq)u^Aq) i^<^.M - n,,(9)) 

{q)^<f>3<P3{q) 

+Vatl (9) (n..3 ('z)n0303 (9) - n,, (q)n,^3 (g)) 
('7)nr,03 (9) - n<,<,3 {q)n,,,,{q)) , 

For the meson loop correction to the quark condensate split (j48p . we have 
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ij 

(ni°j(<7) - IV^Zm) Tt [S+{p)S+{p)S+{p + q)- S-{p)S-{p)S-{p + q)]„^,,^^ + 

2Uaa^{q) Tr [S+{p)S+{p)S+{p + q) + S-{p)S^{p)S-{p + (7)]^„^o - 

J2 nt"^(g) Tr[S+ip)S+ip)S^ip + q)-S^ip)S^{p)S+{p + q)]„^Sm)- 

- n'^Zs^q)) Tr [S+{p)S+{p)S+{p + q) - 5_ (p + q)]^^,,^^ - 
2n,03(g) Tr [S+{p)S+{p)S+{p + q) + + + 

^ n'^^^iq) Tr[S+{p)S+{p)S^{p + q)-S^{p)S^{p)S+{p + q)]^^,^^ 

} 

In complete analogy we may evaluate the diagrams shown in the Figure ^ and get the following explicit vertices 



a. Diagram #3 

Left part: 



i(3) 



E 



E 

1^=1,2 









-q)f{p + k)f{p- 




-k) X 


{ 
















Tr [S+{p)S+{p + k)S+{p- 


^q + k)S+{p + q)- 


f 5- 














^9)] + 






Tr [S+ip)S+ip + k)S+ip- 


^q + k)S+{p + q) - 


-S- 




-(PH 
















Tr [S+{p)S.{p + k)S^{p- 


^q + k)S+{p + q)- 


f 5_ 




hfc)5+(pH 










(n,,,,(fc) - n0303 


(fc)) 


Tr [S+{p)S+{p + k)S+{p- 


^q + k)S+ip + q)- 


f 5- 












-(pH 




2n^03(fc) 




Tr [S+ip)S+{p + k)S+{p- 


^q + k)S+{p + q) - 


-5- 


-{P)S- 








-k)S- 




^q)] - 






Tr [S+{p)S-{p + k)S-{p- 


^q + k)S+ip + q)- 


f 5- 










-k)S- 







E 

i±=l,2 



E 









-q)f{p + k)f{p- 


f 94 


-k) X 


{ 












(ncr£r(fc) no-3(j3 




Tr [S+{p)S+{p + k)S+{p- 


^q + k)S+ip + q) 




-{p)S^ 














2n<,<,3(fc) 




Tr [S+ip)S+ip + k)S+ip- 


^q + k)S+ip + q) 


+ 5- 


^ip)S- 




-fc)5_(pH 










n..(fc) 




Tr [S+{p)S-{p + k)S-{p- 


^q + k)S+ip + q) 


-S- 


-ip)s, 


-(PH 




hgH 




-(PH 




(nr,r,(fc) - n0303 


(k)) 


Tr [S+{p)S+{p + k)S+{p- 


^q + k)S+{p + q) 


-S- 


-{P)S- 


-(PH 










^q)] 


2n^03(fc) 




Tr [S+ip)S+ip + k)S+ip- 


^q + k)S+ip + q) 


+ S- 


-{p)S- 








-k)S- 






n0,(fc) 




Tr [S+ip)S^ip + k)S^ip- 


^q + k)S+ip + q) 


-5_ 


-{p)s, 


4ph 


hfc)5+(pH 




-k)S- 







} 

Right part has been evaluated in (|54l55p . 



b. Diagram #4 



Left part has been evaluated in (|52l53p . 

Right part: 



c. Diagram #5 

Left part: 



(11(7(7 (^) 11(73(72 




~k)S+{p)S+{p^ 


-q)- 


f 5- 


\p)S-{p + k)S- 






^q)] 




Tr[S+ip)S+ip^ 


-k)S+{p)S+{p^ 




-5- 


\p)SAp + k)S- 


-{P)S- 








Tr [S+{p)S-{p^ 


~k)S+{p)S+{p^ 






\p)S+{p + k)S- 


-ip)S- 


(PH 


^q)] 




Xk)) Tt[S+{p)S+{p^ 


^k)S+{p)S+{p^ 




f 5- 


\p)S-{p + k)S. 


Xp)S- 


-(PH 


^q)] 




Tr [S+{p)S+{p^ 


- k)S+{p)S+{p^ 


-q)- 


-5- 


\p)S-{p + k)S^ 


-{p)S- 


-(PH 


-q)] 




Tr[S+{p)S-{p^ 


-k)S+{p)S+{p^ 


^q)- 


f 5- 


^{p)S+{p + k)S- 


-ip)S- 


-(PH 


~q)] 



Lf{q) = ^8nJ ^M'f{p)f{p + k)f{p + q) X { 



E 

i± = l,2 





(A:)) Tr[S+{p)S+{p^ 


-k)S+{p)S+{p^ 


^q) 


-S- 


-(P)^-(PH 


-k)S- 


-ip)S- 


(PH 


-9)] 




Tr [S+ip)S+ip^ 


-k)S+{p)S+ip^ 


^q) 


+ S- 


-(p)5_(pH 




-ip)S^ 


-(pH 






Tr [S+{p)S.ip^ 


-k)S+{p)S+{p^ 


^q) 


-S- 


-(P)^+(PH 


- 


ip)S^ 


(pH 




(n^^(fc) - n0303 


(fc)) Tr[S+{p)S+{p^ 


- k)S+{p)S+{p^ 


'-q) 


-S- 


-(p)^-(pH 


'-k)S- 


-{p)S- 


-(pH 




2n^03(/c) 


Tr[S+ip)S+ip^ 


-k)S+ip)S+ip^ 


-q) 


+ S- 


-(p)5-(pH 


-k)S- 


-{p)S- 


-(pH 






Tr [S+{p)S-ip^ 


-k)S+{p)S+ip^ 


^q) 


-S- 


-(P)5+(PH 


^ k)S- 


-ip)S- 


-(pH 





E 

1^ = 1,2 

} 

Riglit part has been evaluated in (|54l55p . 



Left part has been evaluated in (|52l53p . 
Right part: 



d. Diagram #6 
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